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Abstract. We show that in periodically driven systems, along with the delta-peak at 
the driving frequency, the spectral density of fluctuations displays extra features. These 
can be peaks or dips with height quadratic in the driving amplitude, for weak driving. 
For systems where inertial effects can be disregarded, the peaks/dips are generally 
located at zero frequency and at the driving frequency. The shape and intensity of the 
spectra very sensitively depend on the parameters of the system dynamics. To illustrate 
this sensitivity and the generality of the effect, we study three types of systems: an 
overdamped Brownian particle (e.g., an optically trapped particle), a two-state system 
that switches between the states at random, and a noisy threshold detector. The 
analytical results are in excellent agreement with numerical simulations. 


1. Introduction 

Fluctuation spectra and spectra of response to periodic driving are major tools of 
characterizing physical systems. The spectra are conventionally used to hnd system 
frequencies and relaxation rates and to characterize fluctuations in the system. For 
example, optical absorption spectra give the transition frequencies of atomic systems 
and the lifetimes of the excited states, and the spectrum of spontaneous radiation is a 
well-known example of the fluctuation (power) spectrum [1]. In macroscopic systems 
the spectra are often complicated by the effects of inhomogeneous broadening. Recent 
progress in nanoscience has made it possible to study the spectra of individual dynamical 
systems. A well-known example is provided by optically trapped Brownian particles 
and biomolecules Els], where the power spectra are a major tool for characterizing the 
motion in the trap 015]. Spectra of various types of individually accessible mesoscopic 
systems are studied nowadays in optics 0 0, nanomechanics and circuit quantum 
electrodynamics, cf. jH], biophysics, cf. pun], and many other areas; the technique based 
on spectral measurements has found various applications, photonic force microscopy 
being a recent example, see Ref. m- 

A familiar effect of weak periodic driving is forced vibrations of the system. When 
ensemble-averaged, they are also periodic and occur at the driving frequency up- They 
lead to a (5-shape peak at frequency cuf in the system power spectrum. However, the 
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driving also modifies the power spectrum away from up. A textbook example is inelastic 
light scattering and resonance fluorescence. In the both cases, the system driven by a 
periodic electromagnetic held emits radiation at frequencies that differ from the driving 
frequency [12]. This radiation is one of the major sources of information about the 
system in optical experiments. 

In this paper we study the spectra of periodically driven nonlinear systems. We 
show that, in the presence of noise, along with the (5-shape peak at the driving frequency 
Up, these spectra display a characteristic structure. We are interested in the regime of 
relatively weak driving, where the driving-induced change of the power spectrum is 
quadratic in the amplitude of the driving, as in inelastic light scattering. 

In view of the signihcant interest in the power spectra of systems optically trapped 
in fluids, we consider systems where inertial effects play no role. In the absence of 
driving the power spectra of such systems usually have a peak at zero frequency. In 
particular it is this peak that is used to characterize the dynamics of optically trapped 
particles. 

For a linear system, like a Brownian particle in a harmonic trap, the (5-shape peak 
at Wi? is the only effect of the driving on the power spectrum. This is because motion 
of such a system is a linear superposition of forced vibrations a.t up and fluctuations 
in the absence of driving. The amplitude and phase of the forced vibrations depend 
on the parameters of the system and determine the standard linear susceptibility 0. 
In nonlinear systems forced vibrations become random, because the parameters of the 
system are fluctuating. The power spectrum of such random vibrations is no longer just 
a (5-shape peak (although the (5-shape peak is necessarily present). The driving-induced 
spectral features away from up result from mixing of fluctuations and forced vibrations 
in a nonlinear system. We note the close relation of these features to inelastic light 
scattering, see [Appendix A 


1.1. Qualitative picture 

The idea of the driving-induced change of the power spectrum can be gained by looking 
at a Brownian particle fluctuating in a confining potential, a typical situation for 
optical trapping. The motion of the particle, after proper rescaling of time and particle 
coordinate g, is described by the Langevin equation m 

g = -U'(g) + /(i). U'(q) = dU/dq, (1) 

where U (g) is the scaled potential and f(t) is thermal noise. If potential U (g) is parabolic 
and the system is additionally driven by a force F cos upt, forced vibrations are described 
by the textbook expression 

(?(^)) = exp{-iupt) + C.C., x(a;) = [f/"(geq) - (2) 

where geq is the equilibrium position [the minimum of 17(g)] and x{^) is the susceptibility. 

For a nonlinear system the potential 17(g) is nonparabolic. Because of thermal 
fluctuations, the local curvature of the potential 17" (g) is fluctuating. Intuitively, one 
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can think of the effect of thermal flnctuations on forced vibrations as if t/"(geq) in Eq. ([2]) 
for the susceptibility were replaced by a fluctuating curvature, see Fig. [1] If the driving 
frequency up largely exceeds the reciprocal correlation time of the fluctuations t~^, the 
fluctuations would lead to the onset of a structure in the power spectrum near frequency 
Up with typical width tf^. The quantity t~^ also gives the typical width of the peak 
in the power spectrum at zero frequency in the absence of driving [for a linear system, 
= f/"(geq)]. 


U(q) 


q 

Figure 1. Sketch of a potential of a nonlinear system near the potential minimum. 
Because of the interplay of nonlinearity and fluctuations, the curvature of the potential 
fluctuates. These fluctuations are shown as the smearing of the solid line, which 
represents the potential in the absence of fluctuations. 



Another effect of the interplay of driving, nonlinearity, and fluctuations can 
be understood by noticing that the periodic force causes a periodic change in the 
system coordinate. For a nonlinear system, roughly speaking, this leads to a periodic 
modulation of the local curvature, and thus of Since t~^ determines the shape of 
the zero-frequency peak in the power spectrum, such modulation causes a change of this 
peak proportional to to the lowest order in F. 

Even from the above simplistic description it is clear that the driving-induced 
change of the spectrum is sensitive to the parameters of the system and the noise 
and to the nonlinearity mechanisms. Explicit examples given below demonstrate this 
sensitivity and suggest that the effects we discuss can be used for characterizing a system 
beyond the conventional linear analysis. After formulating how the power spectrum 
can be evaluated in Section |2l we demonstrate the effects of the interplay of driving 
and fluctuations for three very different types of nonlinear systems: an overdamped 
Brownian particle (Section |3]), a system that switches at random between coexisting 
stable states (Section 0]), and a threshold detector (Section [5]). All these systems are of 
broad interest, and all of them display a signihcant driving-induced change of the power 
spectrum. 
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We consider fluctuating systems driven by a periodic force F cos(upt) and assume that 
fluctuations are induced by a stationary noise, like in the case of an optically trapped 
Brownian particle, for example. After a transient time such system reaches a stationary 
state. The stationary probability distribution of the system with respect to its dynamical 
variable q, is periodic in time t with the driving period Tp = 27i/up. The two- 

time correlation function {q{ti)q{t 2 )) [(•) implies ensemble averaging] is a function of 
ti — t 2 and a periodic function of t 2 with period Tp. The power spectrum usually 
measured in experiment is of the form 



(3) 


The correlation function in Eq. (|3]) can be expressed in terms of Pst{q,t) and the 
transition probability density p{qi,ti\q 2 ,t 2 ) that the system that was at position q 2 at 
time t 2 is at qi at time ti >t 2 , 



time-independent coefficients, whereas p{qi,ti\q 2 ,t 2 ) can be expanded in Fexp{±iojpt 2 ) 
with coefficients that depend on ti — t 2 - Therefore the power spectrum ((31) does not 
have terms linear in F. To the second order in F for ca > 0 we have 


$(a;) — $o(^) + — ujp) -|- F'^^p(u!). 


(5) 


The term describes the power spectrum of the system in the absence of driving. 


The term oc 6{u} — up) describes the conventional linear response, cf. Eq. ((2]). However, 
the expression for the susceptibility x{^) nonlinear systems is far more complicated 
than Eq. ([2]); generally, the susceptibility is determined by the linear in F term in 
Pst{q,t). In the optical language, the term oc S{u — up) in ([5]) corresponds to elastic 
scattering of the field F cos upt by the system. 

Of primary interest to us is the term $^ 7 ( 0 ;). This term is often disregarded in 
the analysis of the power spectra of driven systems, while the major emphasis is placed 
on the (5-function in Eq. ([5]). Function <hir(a;) describes the interplay of fluctuations 
and driving in a nonlinear system beyond the trivial linear response. In the considered 
lowest-order approximation in the driving amplitude, does not contain a (5-peak at 
2ujp. However, it may contain a (5-peak at w = 0, which corresponds to the static 
driving-induced shift of the average position of the system. In what follows we do not 
consider this peak, as the static equilibrium position can be measured independently. 

Function ^p can be found from Eq. (jl]) by calculating the transition probability 
density and the stationary probability distribution. This can be done for Markov 
systems numerically and also, in the case of weak noise, analytically, see Secs. [3] and 
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m Alternatively, function can be related to fluctuations of linear and nonlinear 
response of the system and expressed in terms of the fluctuating linear and nonlinear 


susceptibility, see Appendix A We emphasize that the nonlinear response has to be 
taken into account when fluctuations are considered even though we are not interested 
in the behavior of the power spectrum near 2uf or higher overtones or subharmonics of 

Up. 


3. Power spectrum of a driven Brownian particle 

A simple example of a system where d)i?(a;) displays a nontrivial behavior is a 
periodically driven overdamped Brownian particle in a nonlinear conhning potential 
U{q), see Eq. ([T]). This model immediately relates to many experiments on optically 
trapped particles and molecules. We will assume that thermal noise f(t) is white and 
Gaussian and that it is not strong, so that it suffices to keep the lowest-order nontrivial 
terms in the potential. 


1 .2 , 1 1 


U{q) = ^nq^ + + • • •, 


{f{t)f{t')) = 2D5{t-t'), (6) 


where D oc kpT is the noise intensity. In the absence of driving the stationary 
probability distribution is of the Boltzmann form, p® oc exp[—t/(g)/T)]. 

For small D and weak driving force equation of motion q = —U'{q) + f{t)+F cos upt 
can be solved directly by perturbation theory in the noise f{t) and in F, as indicated 


m 


Appendix A Here we develop a different method, which is particularly convenient if 


one wants to go to high orders of the perturbation theory in D and F. 

3.1. Method of Moments 

Systems in which fluctuations are induced by white noise can be studied using the 
Fokker-Planck equation 

dtp = -dq{[-U'{q) FFcosupt] p} + Ddlp. (7) 

This equation can be solved numerically. A convenient analytical approach is based on 
the method of moments, which are dehned as 




dqq"^ / dq'p{q,t + t'\q',t')q'pst{q',t' 


Mn{uj]t') = / dte''' 

Jo J J 

From Eq. ([5|), the power spectrum is <I)(a;) = (2/ri?)Re dt'Mi{u);t'). 

The moments Mn satisfy a set of simple linear algebraic equations 


( 8 ) 


— iuMn{u) + nlF[Mn{uj)] = Dn{n — l)M„_2(a;) 

+ — F nMn-i{u} + iOp) + e nMn-i{u} — up) +Qn+i{t'). (9) 

Here, we skipped the argument t' in and introduced function F[M„] = nMn + 
jJMn+i + ldMn+ 2 -- Functions 


Qn{t) = / 


( 10 ) 
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in the right-hand side of Eq. dH]) can themselves be found from a set of linear equations 
similar to ([9]). They follow from Eq. ([7]), if one sets p = and takes into account 

that pst{q,t) is periodic in t. To the lowest order in F it suffices to keep in Qn(t) only 
terms that are independent of t or oscillate as exp(±ia;irt); respectively, in Eq. (ITOD 


Qn(t) ^ Qn^ + Qn^ exp^iojpt) -f C.C. 


and 


= D{n - l)Qf 2 + FReQ<P„ 

ihipQ'p + nT[QP] = Dn{n - 1 ) 0™2 + ( 11 ) 

The system of coupled linear equations for the moments Mn and Qn can be quickly 
solved with conventional software to a high order in the noise intensity D. Nontrivial 
results emerge already if we keep terms oc DF^\ these are the terms that contribute to 
the power spectrum <hp’(a;) to the lowest order m D. To hnd them it suffices to consider 
terms with n < 3 and Qn with n < 4. This gives 

^ ^ \ ^ [k^ + {u -UFy][K‘^ + {u + upy] 

-374 • ( 12 ) 


This expression refers to |a;| > 0; function contains also a (5-peak at a; = 0, which 

comes from the driving-induced shift of the average static value of the coordinate. 

The solution of the equations for the moments in the considered approximation gives 
a correction oc to the power spectrum in the absence of driving $ 0 ( 1 ^) • To the lowest 
order in D this function displays a Lorentzian peak at a; = 0, = ‘2D/{K^+uy. This 

peak is used in the analysis of optical traps for Brownian particles HE]; with account 
taken of the term oc the zero-frequency peak of ‘ho(<^) becomes non-Lorentzian. 


3.2. Power spectrum for comparatively large driving frequency 

The interpretation of Eq. (IT^ is simplihed in the case where the driving frequency 
exceeds the decay rate, ojp S> n. In this case, periodic driving leads to two well-resolved 
features in the spectrum One is located at w = 0 and has the form 

(2i7/ci;|)(2/3^ — (13) 

This equation can be easily obtained directly by solving the equation of motion 
q = —U'{q)+F cosupt+fit) by perturbation theory in which q{t) is separated into a part 
oscillating at high frequency ojp (and its overtones) and a slowly varying part. To the 
lowest order in F and D, the fast oscillating part renormalizes the decay rate of the slowly 
varying part of q(t), with k ^ k — {F/upy — ( 3 / 2 ) 7 )]. Using this correction in 

the expression for the power spectrum of a linear system $Q°^(a;) = 2D/{nf -|- cu^), one 
immediately obtains Eq. (IT^ to the leading order in n/ujp. 

Interestingly, Eq. flT3|) describes a peak or a dip depending on the sign of 2{3^ — 3'jk,. 
That is, the sign of ^p is determined by the competition of the cubic and quartic 
nonlinearity of the potential U{q). This shows high sensitivity of the spectrum to the 
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system parameters. The typical width of the peak/dip of near a; = 0 is k; the shape 
of the peak/dip is non-Lorentzian. 



CJiK 

Figure 2. Scaled driving induced terms in the power spectrum of an overdamped 
Brownian particle moving in the quartic potential U{q) given by Eq. ([5]), = 

10'^k^^p{u!)/2D. Panels (a), (b), and (c) refer to the scaled cubic nonlinearity 
P'^D/tc' = 0.002 and quartic nonlinearity ^D/iF = 0.0006, 0.00147, and 0.002, 
respectively. The black dots and red solid curves correspond to the numerical 
simulations and Eq. The scaled driving frequency is ujfIk = 5 and the driving 

strength is nF'^/ujpD = 20. For this driving strength and the noise intensity, the 
simulation results in panels (b) and (c) deviate from the theoretical curve. The 
deviation decreases for weaker driving. This is seen from the simulation data in 
panel (b) that refer to kF'^/ uxfD = 5 (blue triangles) and 1.25 (green squares). The 
corresponding spectra are scaled up by factors 4 and 16, respectively. 


The other spectral feature is located aX up and near the maximum has the form 
of a Lorentzian peak, <hjr(a;) ~ {D/ujp)[nf + (cn — ujpY]~^. The height of this peak 
is smaller by a factor jurf 1 than the height of the feature near cu = 0. We note 
that the height of the peak at cjir is proportional to the squared parameter of the cubic 
nonlinearity of the potential /5, but is independent of the quartic-nonlinearity parameter 
7 , to the leading order in the noise intensity D. 

In Fig. [2] we compare the analytical expression (IT^ with the results of numerical 
simulations. The simulations were done by integrating the stochastic differential 
equation q = —f/'(g) + /(t) + Fcosa;i?t using the Heun scheme (cf. [IS]). Panel (a) shows 
that the cubic nonlinearity of the potential leads to a peak at a; = 0 and a comparatively 
small peak aX up. The spectrum becomes more interesting in the generic case where 
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both cubic and quartic terms in the potential are present and is comparable to 7 K. 
Here, as seen from panel (b), as a result of the competition between these terms, 
can have a dip at cn = 0 and two peaks, one near ujp and the other with the position 
determined by and ujp/n. Where the quartic nonlinearity dominates, 7 K S> /9^, 

see panel (c), it is hard to detect the peak aX ujf for small noise intensity. Our analytical 
calculations and numerical simulations show that, for larger noise intensity, this peak 
becomes more pronounced. 

A signihcant deviation of simulations and the asymptotic expression (1121) in panel 
(b) for small a; is a consequence of the near compensation of the contributions to <hir(ci;) 
from the cubic and quartic nonlinear terms in U{q) to the lowest order in and D. 
The terms of higher-order in D and become then substantial. Panel (b) illustrates 
how the difference is reduced if F^ is reduced. We checked that by reducing also the 
noise intensity we obtain a quantitative agreement of simulations with Eq. (1T2|) . 

In some cases the conhning potential of an overdamped system has inversion 
symmetry, and then = 0 in Eq. ([ 6 ]). In such cases spectral features of $ir at the 
driving frequency are oc (yF)^. They can be found by solving the equations for the 
moments with n < 5 and Qn with n < 6 or by solving the equations of motion by 


perturbation theory to the second order in 7 , cf. Appendix A 


4. Power spectrum of a driven two-state system 

We now consider the effect of driving on a two-state system. Various types of such 
systems are studied in physics, from spin- 1/2 systems to two-level systems in disordered 
solids to classical Brownian particles mostly localized at the minima of double-well 
potentials. We will assume that the system dynamics are characterized by the rates 
Wij of interstate i ^ j switching, where = 1,2. In the case of quantum systems, 
this means that the decoherence rate largely exceeds in other words, the typical 
duration of an interstate transition is small compared to \/Wij. For classical systems, 
this description means that small fluctuations about the stable states are disregarded. 


4-.1. The model: modulated switching rates 

A major effect of periodic driving is modulation of the switching rates. It can be quite 
strong already for comparatively weak driving. Indeed, if the rates are determined by 
the interstate tunneling, since the field changes the tunneling barrier, its effect can be 
exponentially strong. Similarly, it may be exponentially strong in the classical limit if 
the switching is due to thermally activated overbarrier transitions, because the driving 
changes the barriers heights. Nevertheless, for weak sinusoidal driving F cos upt the 
modulated rates W}p{t) can still be expanded in the driving amplitude, 

W}p ^ Wij - aijF cosuipt, i,j = 1, 2. (14) 

This equation is written in the adiabatic limit, where the driving frequency up is small 
compared to the reciprocal characteristic dynamical times, like the imaginary time of 
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motion under the barrier in the case of tunneling [12] or the periods and relaxation 
times of vibrations about the potential minima in the case of activated transitions. The 
rates Wij are also assumed to be small compared to the reciprocal dynamical times. 
The driving frequency Wi? is of the order of W^j. 

Parameters aij in Eq. flTTl) describe the response of the switching rates to the 
driving. They contain factors ~ Wtj. Indeed, for activated processes Wtj oc 
exp{—AUi/kBT), where Af/j is the height of the potential barrier for switching from the 
state i. If F cos upt is the force that drives the system, then a^j ~ Wijdi/kpT, where di 
is the position of the ith potential well counted off from the position of the barrier top 
DU- The terms oc which have been disregarded in Eq. flTTll . are oc Wij{di/kBTY in 
this case; a part of these terms that are oc cos2upt do not contribute to $ir(a;) to the 
second order in F, whereas the contribution of the time-independent terms oc comes 
to renormalization of the parameters Wij in $ 0 ( 1 ^); see below. For incoherent interstate 
quantum tunneling, aij oc W^j, too. 

We will use quantum notations |i) {i = 1,2) for the states of the system. One 
can associate these states with the states of a spin-1/2 particle by setting |1) = | )') 
and 12) = | J,). The system dynamics is most conveniently described by the dynamical 
variable q dehned as 

?=|1)(1|-|2)(2| = <t., (15) 

where is the Pauli matrix. For a particle in a double-well potential, q is the coordinate 
that takes on discrete values 1 and —1 at the potential minima 1 and 2, respectively. 

The power spectra of driven two-state systems have been attracting much interest 
in the context of stochastic resonance, see [niiisiisoiini for reviews. By now it has 
been generally accepted that, for weak driving, the power spectrum of a system has a 5- 
peak at the driving frequency with area oc F^, which is described by the standard linear 
response theory [22]. This peak is of central interest for signal processing. However, as 
we show in this Section, along with this peak, the spectrum has a characteristic extra 
structure, which is also oc F^, to the leading order in F. 

4-2. Kinetic equation and its general solution 

It is convenient to write the analog of Eq. (|1|) for the correlation function of the discrete 
variable q as 

='^{iWzp{tM(TzPst{t2)\j) (16) 

Here, p(ti\t 2 ) is the transition density matrix, p(ti|^ 2 ) = X] 1^2)01, and Pst = 

K)(Pst)M(*| is the stationary density matrix. By construction (in particular, because 
of the decoherence in the quantum case) the stationary density matrix is diagonal. Its 
matrix elements (pst)ii give the populations of the corresponding states and periodically 
depend on time, pst(t + 2%/ojp) = Pstif). The transition matrix elements Pij(ti\t 2 ) give 
the probability to be in state i at time ti given that the system was in state j at time 
t 2 . At equal times we have p{t 2 \t 2 ) = I, where / is the unit matrix. 
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Equation does not have the form of a trace over the states |i); rather it expresses 
the correlator in terms of the joint probability density to be in state \j) at time ^2 and 
in state \i) at time ti, with summation over i,j [23]. In the quantum formulation, 
the applicability of this expression is a consequence of the decoherence and Markovian 
kinetics. 

Matrix elements pij{t\t') satisfy a simple balance equation, which in the presence 
of driving reads 

dtpij{t\t') = -WiP{t)pij + wP\t)p 2 j, pij + p 2 j = 1, (17) 


where j = 1,2. Equation for the matrix elements of Pst{t) has the same form, except 
that subscript j has to be set equal to the first subscript. 

From Eqs. flTHl) and flTTll we obtain a general expression for the correlator of interest. 


{(liti)qit 2 )) = exp 
rh 


rti 


dtwp\t) 


X exp 


Jt 2 

dt'wp\t') 


+ {o'z{t2))st I dt\Wp\t) 


't 2 


wr(t)^w:,r(t)±ws\t). ( 18 ) 


r{F) 


F), 


Here, (cr 2 (t))st = (?(i))st = Tr [a^pstit)] is the time-dependent difference of the state 
populations in the stationary state. Generally , {o'z{t))st is nonzero even in the absence 
of driving unless the switching rates are equal, W 12 = H^i. In the presence of driving 
there emerges a periodic term in (cr^(t))st, which describes the linear response, for weak 
driving. 

Disregarding terms oscillating as exp(±2ia;i?f), to the second order in F we obtain 
from the balance equation flTTl) written for (pst)^ 


W P n 

{(Tz{t))st ~ ^ ^ + C.C.] + ^^Rexi(a;F), 

Xi(a;) = 2 (Q!i2hh2i — ot2iWi2) / [hE+ (hF+ — iuj)]. (19) 


Here we introduced notations 


a± = 021 ± hE± = ^21 ± 1Ei2- (20) 

Function Xi(a;) gives the linear susceptibility. In the case of thermally activated 
transitions, Eq. flT^ for Xi coincides with the classical result [1^. The term hE_/lT+ 
gives the difference of the state populations in the absence of driving, whereas the term 
(X gives the time independent part of the driving-induced correction to this difference. 


4 . 3 . The driving-induced part of the power spectrum 


Equation flTSll allows one to calculate the period-averaged correlator {{q{ti)q{t 2 ))) in the 
explicit form and to obtain the power spectrum. As before, we will not consider the 
(5-peak in $( 0 ;) for a; = 0. The spectrum is an even function of a;, and we will consider 
it for a; > 0: 


$ 0 ( 1 ^) 


^H^i2fE2i w+ 

Wf Wl+u^' 
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$ 


^"\uj) = a+ ^ (j)F{iauj,iyujF), 


/i,i/=± 


4>f{.^i ^f) = ~\W+ ~ ~ 


-1 


Q!+11^12W^21 

ulWl 


+ i: 


W_ 




( 21 ) 


2 uf^^+ 

The term $o is the familiar power spectrum of a two-state system in the absence of 
driving HZI. It has a peak at cn = 0 with halfwidth W+ equal to the sum of the 
switching rates. The term describes the driving-induced modihcation of the peak 
centered at a; = 0, 


= K/ 24 )$o(a;) - \xi{ujf)?Wf/{W l + uj^). (22) 

Of major interest to us is the part of the driving-induced term in the power 

spectrum (IMi) . For cn > 0, it shows a resonant peak (or a dip, depending on the 
parameters) at the driving frequency ujf- In contrast to the 5-peak of the linear response, 
the peak has a hnite halfwidth ~ hF_|_ = hFi 2 -l- IF 21 . It is well separated from the peak 
at cj = 0 for ojf S> IF+ and generally is of a non-Lorentzian shape. We stress that, to 
the order of magnitude, the peak has the same overall area as the 5-peak of the linear 
response (in the case of a dip, the absolute value of the area should be considered). 
Another important feature of the peak/dip seen from Eq. fl2T]) is that it is proportional 
to the parameter = 012 + ct 2 i- This parameter describes the change of the sum of 
the switching rates due to the driving. 

For activated switching between potential minima considered in the classical 
stochastic resonance theory, «+ = {kBT)~^{Wi 2 di + W 2 id 2 )- For a symmetric potential 
= 0, since IF 12 = IF 21 and di = —d 2 . Then = 0, in agreement with [21] where a 
symmetric potential was considered. On the other hand, for strong driving it was found 
[25] that the power spectrum for an asymmetric potential displays peaks close to odd 
multiples of the driving frequency and dips close to even multiples of driving frequency. 
In our weak-driving analysis we do not consider peaks/dips near the overtones of uf] 
however, as seen from Eq. CD. the sign of $^^(0;) near uf can be positive or negative, 
depending on the parameters. 

Examples of the driving-induced spectra <I)p’(a;) are shown in Fig.[3l One can clearly 
see the peaks or dips both at ca = 0 and at the driving frequency uf- In agreement with 
Eqs. (j2T|) and (j22j) . the signs of the features of <I>i? are determined by the interrelation 
between the parameters of the two-states system. For illustration purpose we chose the 
values of the ratio of the response parameters a 2 i/Q!i 2 to he between plus and minus the 
ratio of the switching rates in the absence of driving, W 2 \IW\ 2 - As seen from Fig. |3l 
the spectra are very sensitive to the ratio a 2 i/«i 2 - We have seen this sensitivity also for 
different values of IF 2 i/hFi 2 . 

Unexpectedly, a hnite-height spectrum ^^^(ci;) emerges even where the linear 
susceptibility is equal to zero, which happens for Q;i 2 hF 2 i = a! 2 ihFi 2 . This is seen from 
Eq. fl2T|) and also from Fig. [3l The red line with Q! 2 i/q!i 2 = 7/3 refers to this case, 
and the area of 5-peak in the spectrum is zero. As seen from the hgure, numerical 
simulations are in excellent agreement with the analytical expressions. 
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Figure 3. The driving induced terms in the power spectrum of the two-state system 
for the ratio of the switching rates W 21 /W 12 = 7/3. The scaled driving frequency 
and amplitude are u}f/W+ = 5 and FctxijWxi = 1. On the thick solid (red), dot- 
dashed (black), long-dashed (blue), short-dashed (green), and thin solid (purple) lines 
the ratio aixjaxi is 7/3, 7/6, 0 , —7/6, and —IjZ. The vertical line at wf shows the 
position of the (5-peak at wf- The areas of the (5-peaks for different aixjciYi are given 
by the heights of the vertical segments. The heights are counted off from the lines to 
the symbols of the same color, i.e., to the circle, triangle, and open and full square, in 
the order of decreasing a 2 i/ai 2 ; there is no symbol for 0 ( 21 /q!i 2 = 7/3 as there is no 
(5-peak in this case. The inset shows the full spectrum with (red) and without (black) 
driving for aixjoLYi = 7/3. The curves and the dots show the analytical theory and 
the simulations, respectively. 


The structure of the spectrum near a; = 0 will be modihed if one takes into account 
terms oc in the expressions for the switching rates flTTl) . In the considered leading- 
order approximation in F these terms have to be averaged over the driving period and 
are thus independent of time. The correction due to these terms can be immediately 
found from Eq. (12T]) for <ho('^) by expanding $0 to the hrst order in the corresponding 
increments of Wij] this correction is of a non-Lorentzian form. 


5. Threshold detector 

An insight into the dynamical nature of the driving-induced change of the power 
spectrum can be gained from the analysis of the spectrum of a threshold detector. 
Such detectors are broadly used in science and engineering, and their analogs play an 
important role in biosystems. We will employ the simplest model where the output of 
a threshold detector is g = — 1 if the signal at the input is below a threshold value rj, 
whereas g = 1 otherwise, and will consider the case where the input signal is a sum of 
the periodic signal F cos upt and noise f(t), 

(/(() = 29 [F(() +q() - t;| - 1, 


(23) 
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where Q{x) is the Heaviside step function. To avoid singularities related to non¬ 
differentiability of the ©-function, we will model the output by 

q{t) = tanh[A{F{t) + f{t) - r])], A > 1, (24) 

and in the hnal expressions will go to the limit A —)■ cxd. Much work on the interplay 
of noise and driving in threshold detectors has been done in the context of stochastic 
resonance, cf. [261 Wi\ 128] . In these papers of primary interest was the signal to noise 
ratio; the issues we are considering here, i.e., the occurrence of the effective “inelastic 
scattering” and “fluorescence” as a result of interplay of nonlinearity and noise, have 
not been addressed, to the best of our knowledge. 

In the absence of noise, the power spectrum of q{t) is a series of (5-peaks at ojf 
and its overtones (including cu = 0), provided the driving amplitude F > rj, whereas 
for F < ?7 we have q = —1 and the power spectrum is just a (5-peak at cu = 0. On 
the other hand, if F = 0 and f{t) is white noise, in the limit A —)■ cx) in Eq. (]2T|) the 
correlator {q(t)q(t')) = 0 for t 7 ^ t', since the values of q(t) at different instants of time 
are uncorrelated and {q) 0 , whereas 1 . 

The singular behavior of the correlator {q{t)q{t')) in the case of white noise persists 
also in the presence of driving. This is a consequence of the absence of dynamics, i.e., 
any memory effects in the variable q{t) ( 12 ^ . and the singular distribution of white noise, 
where the intensity (^^(t)) diverges. 

Dynamics can be brought into the system by the noise color. Such noise can be 
thought of as coming from a dynamical system with retarded response, which is driven 
by white noise. We will be interested in the correlator {q{t)q{t')) and the power spectrum 
$( 0 ;) for weak driving, where the driving amplitude is F 77 (subthreshold driving), 
and for a simple colored noise, the Ornstein-Uhlenbeck noise. This is Gaussian noise 
with correlator 


{at)i(t')) = {D/K)exp(-K\t-t'\). 


(25) 


Parameter k characterizes the decay rate of noise correlations. 

Because the threshold detector has no dynamics on its own, the value of the 
variable q{t) is determined by the instantaneous value of the noise ^{t). We can write 
q{t) = q{t,f{t)), where q{t,f) is given by Eqs. fl23ll . fl2Tl) with f^{t) replaced with C,. Then 
the general expression for the correlator of q{t), Eq. (jl]), can be rewritten as 

(9(G)g(G)) = j dfidf2q{tl,^l)q{t2,^2) (26) 


Here, the superscript f indicates that the corresponding transition probability density 
and the stationary distribution refer to the random process ^{t). 

The form of the transition probability for the process (l25l) is well-known [29], 

k(6 - 




K 


27rF(l — 


exp 


2F(1 - 


(27) 


The stationary distribution PsP(^iGi) is given by the same expression with ^2 —t — 00 . 
Substituting these expressions into Eq. (12B]) and expanding q(t,f) in F(t), after averaging 
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over the driving period we obtain to second order in F{t) for ti > t 2 


{{(l{ti)q{t2))) = C + 4: I dfi I df2 -pSP(6) PsP(6) 

J T) J T) ^ 

+ 2F‘^ cosujpiti - t2)p^^\r],ti\7],t2)pfX'n) 


SO, 


- 2F^ 


df2pTi^2)^ P^^Hv,tl\^2,t2) - p’i\v) 

dp L 


(28) 


Here, C is a constant independent of time; it leads to a <5 peak at a; = 0 in the power 
spectrum and will not be considered in what follows. The remaining terms are time- 
dependent. They decay with increasing |fi — ^ 2 !, except for the term that oscillates as 
exp[±za;ir(fi — ^ 2 )] and describes the standard linear response to periodic driving. As 
seen from Eq. fl25]) . this term has the form 

2F^COSUJF{ti-t2) p^iv) =^F^\x{ujf)\'^ COSCJpiti -t2), 

x{u) = 2pif{p) ^ {2n/7rDY/^exp[-Kpy2D), (29) 

where x(a;) is the standard linear susceptibility [I3] of the threshold detector. 
Interestingly, this susceptibility is independent of frequency. This is because the detector 
has no dynamics, its response to the driving is instantaneous. An alternative derivation 
of the expression for the susceptibility, which provides a useful insight into the response 
of the threshold detector, is given in Appendix A It also shows how to deal with the 
singularities in Eq. fl28|) for ti —)■ ^ 2 , which emerge after the transition A —)■ cx) in 
Eqs. fl2T|) and fl26|) . 

The power spectrum $( 0 ;) is obtained from Eq. (|28ll by a Fourier transform. The F- 
independent term in Eq. (|28|) gives the power spectrum $ 0 ( 1 ^) in fhe absence of driving. 
It has a peak at cj = 0. The term cx cosa;ir(fi — ^ 2 ) gives a h-peak and also a finite-width 
peak F‘^^p\uj) at frequency ojf- The last term in Eq. fl28l) gives a driving-induced 
feature in the power spectrum at zero frequency 

The shape of the spectra is determined by the dimensionless parameter that 
characterizes the ratio of the threshold to the noise amplitude p{k/DY^^. For weak 
noise, where ^ 1, the peak near ujf has the form 

- 1/2 




1 




Re 


np 


+1 


.UJ — OJp 


^—KTf I2D 


(30) 


AD K, 

Here we assumed that wf/k is sufficiently large, so that the features of d’i? centered at 
ujf and a; = 0 are well separated; Eq. fl5U]) applies for \uj — ujf\ ojf- The spectrum 
fl30|) has a characteristic non-Lorentzian form with typical width n^p^/AD. However, its 
area is small. 

In the opposite limit of low threshold, p{k/DY^'^ 1, to the leading order 


«>«(a.) 


2^D 


Re 


r 


.OJ 


2k 


"A/r 


.UJ 


ujf 


2k 


(31) 


near uf- This spectrum falls off slowly away from the maximum, as \uj — ujf\ for 
Iw — oii?! S> K. Equation flTI]) does not contain the threshold p. The small-r; correction 
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to (1^ for a; = Wi? is (1 —\n 2 )K,r]'^/ 7 iD‘^. It is positive. From the comparison of Eqs. (I5U]) 
and fl3T|) . one sees that the height of the peak at up hrst increases with the increasing 
r]{K/DY/^, but then starts decreasing. 

In Fig. |4]we show analytical results for the power spectra obtained from Eq. fl28|) for 
several parameter values and compare them with the results of simulations. Immediately 
seen from this hgure is that the driving modihes the overall spectrum most signihcantly 
near cu = 0 and near up for large ojp/n. There emerges a hnite-width peak at ojp. As 
seen from the inset in panel (b), the width of this peak increases with decreasing noise 
intensity, that is, with increasing This is a counterintuitive consequence of 

the unusual interplay of noise and driving in a threshold detector. The height of the 
peak displays a nonmonotonic dependence on rj^K/ 



Figure 4. Power spectrum of the threshold detector, (a): The full power spectrum; 
the scaled frequency and the intensity of the driving are = 100 and F^n/D = 

0.0025. The scaled threshold is rjlnjDY^'^ = 0.5. Inset: the spectrum near the driving 
frequency. The delta peak has been subtracted. The curves and black dots refer to 
the theory and simulations, respectively, (b): The low-frequency part of the driving- 
induced term in the power spectrum for wp/n = 50 as given by Eq. (E51) . The 
solid (black), long-dashed (red), short-dashed (blue) and dot-dashed (green) curves 
correspond to the scaled value of the threshold ylnlDY^'^ = 0.1,0.8,1.2, and 2. Inset: 
the spectrum near the driving frequency, uipln = 50. 

The low-frequency spectrum ~ $^^( 0 ;) also displays a pronounced feature 

near a; = 0. One can show from the analysis of the last term in Eq. (1281) that, for 
small ri{n/DY''^, this feature is a dip, with $^^(0) = —1/D for rj^n/DY^'^ —^ 0. The 
shape of the dip is non-Lorentzian, with typical width k. As ti^k/DY^'^ increases, the 
depth of the dip decreases. Ultimately the shape of the spectrum changes completely. 
For large rj^n/DY^^ the spectrum becomes broad and shallow. To the leading order 
in [ri{K/DY^Y~^y t)e written as {2 / ti D){D/exp{—/2 D)^p\2Du/n'^riY , 
where the dimensionless function is zero for a; = 0, has a minimum at x ~ 1.7, 

where it is ~ —0.6, and then approaches zero with increasing x as x~^/‘^. 
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6. Conclusions 

The results of this paper demonstrate that the interplay of driving and fluctuations leads 
to the onset of specihc spectral features in the power spectra of dynamical systems. 
Such features are analogs of inelastic light scattering and fluorescence in optics, where 
an electromagnetic held can excite radiation at a frequency shifted from its frequency 
and also at the characteristic system frequency. Our results show that, in classical 
systems and in incoherent quantum systems, the spectral features emerge as a result of 
the huctuation-induced modulation of the response to the driving. Such modulation is 
common to nonlinear systems. 

Since nonlinearity and noise are always present in real systems, the occurrence 
of the driving-induced spectral features in the power spectra should be also generic. 
However, these features are specihc for particular systems, which allows using them for 
system characterization. 

We have studied three types of systems, all of which are attracting signihcant 
interest in mesoscopic physics and in several other areas of science. The hrst one is an 
overdamped Brownian particle huctuating in a nonparabolic potential well. This model 
describes, in particular, small particles and molecules optically trapped in a liquid. We 
hnd that, when the particle is periodically driven, the nonparabolicity of the potential 
leads to an extra spectral peak or a dip at zero frequency. For comparatively weak noise, 
the sign of the driving-induced term in the spectrum at small u is determined by the 
competition of the cubic and quartic nonlinearity of the potential. The overall shape 
of the low-frequency spectrum strongly depends on the form of the conhning potential 
as well. In addition, along with a (5-peak at the driving frequency, the driving-induced 
spectrum displays a peak at this frequency with a width of the order of the relaxation 
rate of the system. 

We have also studied a two-state system that at random switches between the 
states. We assumed that the driving modulates the rates of interstate switching. The 
driving-induced spectrum has a rich form. Depending on the interrelation between the 
switching rates without driving and the driving-induced corrections to the rates, it can 
have peaks or dips both at cn = 0 and at the driving frequency. The typical width of 
the peaks/dips is given by the sum of the interstate switching rates without driving. 
Interestingly, these hnite-width spectral features can emerge even where the (5-peak at 
the driving frequency has very small (or zero) intensity. 

The third system we studied is a threshold detector. Here the dynamical nature 
of the driving-induced spectral change is particularly pronounced, as this change does 
not occur if the noise in the detector is white, except for the (5-peak at the driving 
frequency. On the other hand, for colored noise driving does change the power spectrum 
nontrivially. As in other systems, we hnd a driving-induced spectral feature near zero 
frequency. It can be a peak or a dip depending on the ratio of the threshold to the 
appropriately scaled noise intensity. Also, the height of the hnite-width peak at the 
driving frequency displays a nonmonotonic dependence on this ratio, as does the width 
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of the peak, too, i.e., noise can both increase or decrease the width. 

In all studied systems inertial effects played no role: the peaks of the power spectra 
are located at zero frequency in the absence of driving. Therefore driving-induced 
spectral features near the driving frequency and zero frequency correspond to inelastic 
scattering and fluorescence, respectively. However, in contrast to the conventional 
fluorescence, driving can induce a dip in the spectrum at zero frequency, as we have 
seen in all studied systems (the total power spectrum remains positive, of course). The 
occurrence of the dip looks as if the driving were decreasing the noise in the system, 
although in fact the dip has dynamical nature. 

The power spectra of weakly damped nonlinear systems should also display extra 
features in the presence of weak periodic driving. The effect should be most pronounced 
where the driving is resonant. Along with the features near the driving frequency and 
near w = 0, there should arise features near the eigenfrequencies of slowly decaying 
vibrations about the stable states. Several features of the power spectra have been 
studied for nonlinear oscillators in the regime of strong driving, see recent papers [32l [33] 
and references therein. Interestingly, the results do not immediately extend to the weak- 
driving regime, and the features of the interplay of nonlinearity and driving where they 
are of comparable strength remain to be explored. However, it is clear from the presented 
results that the driving-induced change of the spectra is a general effect that provides 
a sensitive tool for characterizing fluctuating systems and their parameters. 

The research of YZ and MID was supported in part by the US Army Research Office 
(grant W911NF-12-1-0235), US Defense Advanced Research Agency (grant FA8650-13- 
1-7301), and by TOYOTA Central R&D Labs., Inc. 

Appendix A. Formulation in terms of fluctuating susceptibilities 

The change of the power spectrum induced by the driving can be analyzed in terms 
of the fluctuating linear and nonlinear susceptibility of the system. If the dynamical 
variable that describes the state of the system is q{t), to the second order in the driving 
Fit) we have 

q{t)^qo{t)+f dt'xi{t,t')F{t') + f f dt'dt"x 2 {t,t',t'')F{t')F{t''), (A.l) 

J —oo J J —oo 

where qo{t) is the (random) value of q{t) in the absence of driving. Functions 
Xi and X 2 describe the linear and nonlinear response. We emphasize that these 
functions themselves are random, there is no ensemble averaging in Eq. (lA.lll . This 
equation is merely a consequence of the causality principle. Spatial and temporal 
fluctuations of the linear susceptibility xi are standardly considered in the context 
of light scattering [301 |3T] . However, the analysis of the power spectrum of nonlinear 
systems in the presence of driving requires also taking into account the fluctuating 
nonlinear susceptibility X 2 - 

The linear and nonlinear fluctuating susceptibilities lead to two terms in the driving- 
induced part of the power spectrum dehned by Eq. (|5|), <I)i?(a;) = ^^^^(a;) -|- $^^(a;). 
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Substituting into flA.ip F(t) = F cosojpt exp(et) with e +0, we obtain [31] 

1 roo r pO 

$W(a;) = -Re/ 

^ Jo JJ-oo 

X (xi(i.« + «')lxi(0,i")-(xi(0,i")>l). (A.2) 

and 

POO P pO 

$p^(ci;) = Re / dte^^ // dt'dt" cos[ujF{t'— t")] 

Jo J J—oo 

X [(X2(t, t + t',t + f)qo{Q)) + (go(^)X2(0, t', t”))], (A.3) 

where we assumed (qoit)) = 0. 

The general form of Eqs. flA.2p and flA.3D immediately shows two distinct effects of 
the driving, which are pronounced where the driving frequency ujp largely exceeds the 
typical relaxation rate of the system. The term is a function of the detuning 

of frequency u from the driving frequency up- Therefore one may expect that 
displays features like peaks or dips near ojf- In contrast, should display features 

near the characteristic frequencies of the system. In particular, for overdamped systems 
that we consider here such features occur around a; = 0. We note that may also 
display features at the system eigenfrequency, since the integrand in Eq. (IA.2D depends 
on Up- 

A convenient way to calculate the fluctuating susceptibilities Xi ,2 is based on solving 
dynamical equations of motion of the system. For example, for an overdamped Brownian 
particle described by the Langevin equation q = —U'{q) + f (t)+F cos upt with nonlinear 
potential ([3|), one can proceed by rewriting this equation in the integral form, 

X [F COS uipt'+ f{t')]. (A.4) 

For small / and F, one can then expand the g-dependent exponential in the right-hand 
side and use successive approximations in F and /. The fluctuating susceptibility xi 
is given by linear in F terms, whereas X 2 is given by the terms quadratic in F. The 
advantageous feature of this method is that it is not limited to white noise. However, 
the method becomes impractical if the noise intensity is not weak, and even for weak 
noise it becomes cumbersome if one goes to high-order terms in the noise intensity. 

We have checked that the calculation based on Eq. flA.4D gives the same result for 
the driving-induced part of the power spectrum <f)p’(c(;) as the method of moments. We 
have also found that, in the second order in the noise intensity D, the term 7g^/4 in 
U{q) leads to the onset of a peak in at up. 

Appendix A.l. Fluctuating susceptibility of a threshold detector 

Fluctuating linear susceptibility has a particularly simple form for a threshold detector. 
By linearizing in F{t) expression (j23|) for the output of the detector, we obtain from the 
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definition of the susceptibility flA.ip 

X{t, t') =25{t-t' - 0)5{f{t) - ri), (A.5) 

where rj is the threshold and f{t) is the noise. Zero in 6 {t — t' — O) reflects causality: the 
detector output q(t) is determined by the value of the driving just before the observation 
time; the very 5-function indicates that the effect of the driving is not accumulated over 
time, the response is instantaneous (but causal). 

The standard linear susceptibility x(cu) is given by expression 

poo 

X(u;) = / 

Jo 

From (lA.Sp . x{^) = where pltip) is the stationary probability density of 

the noise f{t), cf. Eq. (EOj). It applies for an arbitrary noise f{t), not just for the 
exponentially correlated noise considered in Sec. [5l 

Similarly, the fluctuating nonlinear susceptibility of the detector is 

X2{t, t', t") = -5{t - t' - 0)5(f - t" - t])dr^6 {f{t) - f]). (A.6) 

Substituting Eqs. (lA.Sp and flA.bp into the general expressions for the power spectrum in 
terms of fluctuating susceptibilities, Eqs. flA.2p and flA.Sp . we obtain the power spectrum 
in the same form as what follows from Eq. fl28|) . 
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